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Abstract 

DFT calculations performed on Si2He, Si2F6, Si2Cl6, and Si2Br6 are reported. 
The evolution of the energy, the chemical potential and the molecular hardness, as a 
function of torsion angle, is studied. Results at the DFT-B3LYP/6-311++G** level 
show that the molecules always favor the stable staggered conformations, with low 
but significant energy barriers that hinder internal rotation. The chemical potential 
and hardness of Si2H6 remains quite constant as the sylil groups rotate around the 
Si-Si axis, whereas the other systems exhibit different degrees of rearrangement of 
the electronic density as a function of the torsion angle. A qualitative analysis of the 
frontier orbitals shows that the effect of torsional motion on electrophilic attack is 
negligible, whereas this internal rotation may generate different specific mechanisms 
for nucleophilic attack. 
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1 Introduction 



Recently disilane has attracted attention (0; 0) due to its importance in the 
production of sihcon based semiconductor devices. Its geometry is quite similar 
to ethane, which is the best known and most widely studied example (0; ^) 
of simple molecules with properties that markedly depend on the rotation of 
a group of atoms around one or more internal bonds, going through stable 
and unstable conformations as a full 360° rotation is executed. In particular, 
the central C-C bond of ethane is a threefold symmetry axis. Thus, as one 
of the two methyl groups rotates around this axis the molecule goes through 
(stable) staggered and (unstable) eclipsed conformations (see Fig. 1). The 
preferred staggered structure is attributed to steric effects (^ ^; ^, 
more precisely to increased repulsion between electrons in bonds that are 
drawn closer together (^. 

On the other hand, the fundamental processes in the disilane decomposition 
on silicon surfaces are relevant to the understanding and optimization of the 
growth of epitaxial silicon on silicon substrates. The morphological parameters 
of the eclipsed and staggered silane conformations were recently calculated by 
Pophristic and co-workers (||). They concluded that the origin of the eclipsed 
to staggered relaxation is related to the preferential hyperconjugative stabi- 
lization (meaning energy stabilization through electron excitation to a delo- 
calized state). This charge delocalization changes the electronic properties of 
the molecule, as a function of the conformation it adopts. 



When a reaction moves forward along the reaction coordinate, a redistribu- 
tion of the ground-state electron density takes place, and the resulting energy 
change can be understood in terms of the response of the system to variations 
of the total number of electrons A^, and of the external f(f^ potential (|ll]). 
Density functional theory (DFT) (|1T]; 0) has been quite successful in provid- 
ing a theoretical basis for qualitative chemical concepts like chemical potential 
(/x) and hardness [r]), which describe the response of the system when is 
varied for a fixed v{r) (pHI). In DFT is the Lagrange multiplier associated 
with the normalization constraint that requires conservation of the number of 
electrons A^. Classical structural chemistry is recovered with the identification 
of /i as minus the electronegativity (// = — x), a well known and well estab- 
lished quantity. Definitions of and rj, two global electronic properties that are 
implied in the reactivity of molecular systems, were given by Parr et al. ([13|) 
and Parr and Pearson (|T1|; ITB|), respectively. The application of DFT concepts 



to the analysis of chemical reactions is better appreciated with the help of the 
principle of maximum hardness (PMH), that asserts that molecular systems 
reach equilibrium tending towards states with the highest hardness (|1^; 
In this context the PMH can also be helpful in identifying transition states 
where minimum values of rj are expected (pISl). 
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The main purpose of this paper is to provide a detailed report on the geo- 
metric and electronic structure of disilane (Si2H6) and the family of closely 
related molecules 812X5, where X = F, CI and Br, as well as the implications 
this structure has on the molecular properties. We focus our attention on the 
changes that are induced on the energy and molecular properties as the molec- 
ular conformation periodically changes from staggered to eclipsed and back to 
staggered through rotation with respect to the Si-Si bond of the SiXs group 
of SiaXg (X=H,F,a,Br). 

This paper is organized as follows: after this introduction we discuss technical 
aspects of our calculation in Sec. 2, discuss the molecular geometry in Sec. 3, 
the electronic energy profiles and rotational barriers in Sec. 4 and the chemical 
potential and the hardness in Sec. 5. In Sec. 6 we present a qualitative analysis 
of the chemical reactivity of silanes and finally, we close the paper in Sec. 7 
drawing conclusions. 




Fig. 1. Staggered (left) and eclipsed (right) conformations. The rotation angle a = 0' 
for the staggered and a = 60° for the eclipsed conformations, respectively. 



2 Computational Details 

The calculations reported in this paper were performed using the Gaussian 
98 (^ package. The results we report were obtained implementing the B3LYP- 
DFT method and corresponds to the B3LYP/6-311++G** level. A lower 
level B3LYP/6-31G* calculation was tested for comparison purposes. For the 
staggered and eclipsed conformations, and in order to test the DFT results, 
MP2/6-311++G** calculations were also carried out to check the quality of 
the B3LYP results. 
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3 Geometry 



Constrained geometrical optimization was performed varying the dihedral an- 
gle a, defined as the rotation angle of the silyl group of Si2H6 (or the three 
X atoms in Si2X6), located at one end of the Si-Si bond, relative to the same 
three atoms at the other end of this bond (see the illustration in Fig. 1). The 
angle < a < 60°, with a — defined as the staggered conformation and 
a — 60° as the echpsed one, is varied in steps of 10°. 

The structural parameters obtained for the staggered conformation are fisted 
and compared to experimental data -whenever the latter is available- in Ta- 
ble 1. It is apparent that for the lighter molecules (i.e. Si2He and Si2Fg) 
B3LYP/6-31G* yields better agreement with experimental values than B3LYP 
/6-311++G**, while the larger basis set 6-311++G** fares better for SisCle. 
We expect the same to hold for Si2Br6 (a molecule that has not yet been syn- 
thesized) since heavier atoms require larger basis sets for a proper description. 
Thus, in what follows below, our comparisons with experiment are based upon 
the results of B3LYP/6-311++G**. 

Rotation of one syhl group with respect to the staggered conformation is 
accompanied by a significant change in the Si-Si distance (see Fig. 2) while the 
Si-X (X=H, F, CI and Br) distance remains almost unchanged (see Table 1). In 
fact, Si2Br6 displays the largest deformation, which amounts to about 1.59%, 
while Si2F6 undergoes a tiny elongation of only 0.23% . 

The angle Z(SiSiX), between the Si axis and the X-atoms, exhibits a small but 
systematic increase as a function of a. Again, this change is largest for Si2Br6 
(approximately 0.66%) and smallest for Si2F6 (approximately 0.04%). 



4 Energy Profiles and Rotational Barriers 

Fig. 3 shows the evolution of the total energy for each molecule studied, mea- 
sured with respect to the total energy in the staggered configuration. In each 
case, staggered conformation is of minimum energy and eclipsed conformation 
presents maximum energy. 

It is also evident from Figs. 2 and 3 that the energy follows the same trend 
as the Si-Si distance along the torsional angle. Clearly, the torsional potential 
energy can be understood in terms of the structural changes of the molecule 
undergoes as a is varied. Si2F6 being almost free to rotate, in the sense that 
it undergoes only minor geometrical changes, presents a rather small rota- 
tional barrier of ~0.61 kcal/mol, while Si2Br6 has a rotational barrier of ~ 2.6 
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Molecule 


B3LYP/6-31G* 


B3LYP/6-311++G** 


Experiment 


Si2H6 (staggered) 








d(Si-Si) 


2.350 


2.354 


2.331 @ 


d(Si-H) 


1.489 


1.487 


1.492 d) 


Z (SiSiH) 


110.4 


110.6 


110.3 d) 


Si2H6 (eclipsed) 








d(Si-Si) 


2.360 


2.366 




d(Si-H) 


1.489 


1.487 




Z (SiSiH) 


110.8 


110.6 




Si2F6 (staggered) 








d(Si-Si) 


2.319 


2.336 


2.317 d) 


d(Si-F) 


1.593 


1.598 


1.564 d) 


Z (SiSiF) 


110.5 


110.7 


110.3 d) 


Si2F6 (eclipsed) 








d(Si-Si) 


2.326 


2.341 




d(Si-F) 


1.592 


1.598 




Z (SiSiF) 


110.7 


110.7 




Si2Cl6 (staggered) 








d(Si-Si) 


2.355 


2.354 


2.320 


d(Si-Cl) 


2.060 


2.056 


2.002 


Z (SiSiCl) 


109.7 


109.6 




Si2Cl6 (eclipsed) 








d(Si-Si) 


2.377 


2.378 




d(Si-Cl) 


2.059 


2.056 




Z (SiSiCI) 


110.0 


109.9 




Si2Br6 (staggered) 








d(Si-Si) 


2.335 


2.368 




d(Si-Br) 


2.211 


2.232 




Z (SiSiBr) 


108.7 


109.2 




Si2Br6 (eclipsed) 








d(Si-Si) 


2.356 


2.405 




d(Si-Br) 


2.209 


2.232 




Z (SiSiBr) 


108.7 


110.0 





Table 1 
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20 40 60 



a(degrees) 



Fig. 2. Change of the Si-Si distance for the four molecules studied, in percentages. 
The open circles are the calculated points and the lines are guides to the eye. 




20 40 60 



a (degrees) 

Fig. 3. Electronic energy as function of torsion angle a. The open circles are the 
calculated points and the lines are guides to the eye. 

kcal/mol, consistent with its larger geometric changes. The above results indi- 
cate that torsional potential barriers, that hinder the interconversion between 
two staggered conformations, arise from structural rearrangements induced by 
an interplay between steric repulsion and hyperconjugation effects (0). 

It is important to remark that at the B3LYP/6-31G* level Si2Br6 is predicted 
to be stable in the eclipsed configuration, with a rather significant energy dif- 
ference of 1.17 kcal/mol relative to the staggered one. MP2/6-31G* also yields 
a smaller energy for the eclipsed configuration, but with a much smaller dif- 
ference of only 0.150 kcal/mol. However, the MP2/6-311++G** calculations 
agree with the B3LYP/6-311++G** results. Thus, it seems that 6-311++G** 
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Molecule 


Rotational Barrier (kcal/mol) 


Si2H6 


0.9441 




0.6096 


SiaCle 


1.6272 


Si2Br6 


2.5920 



Table 2 

Calculated rotational barrier at the B3LYP/6-311++G** level for 
Si2H6,Si2F6,Si2Cl6and Si2Br6 

is the minimum basis set required to correctly describe the rotational behavior 
of Si2Br6. 

Table 2 displays the B3LYP/6-311++G** rotational barriers we obtained. 
Substitution of the hydrogens, by the more electronegative fluor atoms, results 
in a lowering of the potential barrier. The electronic population is now mainly 
localized at the SiFa groups thus weakening the Si-Si torsional bond. The 
values for Si2H6 and Si2F6 compare well with those calculated by Cho et 
al. (^. The experimental values for the rotational barrier of Si2H6 are ~ 1 
kcal/mol, and for Si2F6 between 0.51-0.73 kcal/mol, according to early electron 
diffraction measurements (0). 

On the other hand, substitution of the hydrogens by chlorine and bromine 
atoms tends to keep the electronic population uniformly distributed, and the 
observed increase of the potential barrier seems to be related to steric hin- 
drance between quite voluminous chemical groups. 



5 Chemical Potential and Hardness 



In DFT the chemical potential of a molecule is defined by the derivative of 
the energy with respect to the number of electrons N at constant external 
potential v{r): 



/^=h^ > (1) 




where E is the energy and the number of particles. For a finite system this 
extrapolation takes the form ([TT|) 

^^hyE{N + l)-E{N-l)]. (2) 
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Moreover, following Koopmans' theorem (|T9[), the anion energy E{N + 1) can 
be approximated by E{N + 1) f» E{N) + -E'lUMO' cation energy 

E{N - 1), by E{N - 1) ^ ^(A^) - %0M0' where ^lUMO and %0M0 
are the energies of the Lowest Unoccupied and Highest Occupied Molecular 
Orbital, respectively. Within this approximation 

^ - 2(^LUM0 + %0M0) • (3) 



Another relevant characteristic property we want to probe is the chemical 
hardness rj, defined as 

\ / v(r) 



which can be approximated by a finite difference as follows: 

rj = ^[E{N + 1) + E{N - 1) - 2E{N)] , (5) 



which in terms of the HOMO-LUMO energies reads 

V = ^[^LUMO - %0M0] • (6) 



Nevertheless it is important to notice that the actual changes in the torsional 
energy must include the geometrical changes induced by the removal or addi- 
tion of electrons. This geometrical relaxation may include symmetry changes, 
such as the rotations around the Si-Si axis we study in this paper. Thus, the 
significance of t] and as calculated above, whether with the HOMO-LUMO 
approximation or with the unrelaxed (or constraint relaxed) cation and anion 
energies, is not completely accurate. 

Figs. 4 and 5 display the chemical potential fi and chemical hardness rj, re- 
spectively, as a function of the torsion angle a, in the HOMO-LUMO ap- 
proximation. An appreciably difference of the n values between two reference 
conformations implies that an electronic rearrangement, with some charge 
transfer from the higher towards the lower /i conformation, will take place. By 
inspection of Fig. 4 we observe that the chemical potential versus a profiles 
for all molecules, except Si2He, display a variation of ~ 1 kcal/mol as a 
varies by 60°, always opposite in sign to the relative to the energy variation 
displayed in Fig. 5. This is an indication that torsion implies a rearrangement 
of the electronic density. In contrast, for Si2He, the chemical potential remains 
quite constant over the range > a > 60°, with A/i ~ 0.10 kcal/mol. 



8 



o 
E 

15 
o 



LU 

+ 

—I 
LU 



O 

E 



LU 

+ 



-97.62 

-97.67 

-97.72 

-97.77 f 

-97.82 
-123.0r 



-123.8 
-124.0 



1 1 1 1 1 


1 ^ \ ^ \ 




Si^Br, 



10 20 30 40 50 80 10 

a(degrees) 



30 50 

a(degrees) 



-140.0 
-140.5 
-141.0 
-141.5 
-142.0 
-125.5 
-126.0 
-126.5 
-127.0 



Fig. 4. Chemical potential in the HOMO-LUMO approximation 
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Fig. 5. Chemical hardness in the HOMO-LUMO approximation. 

Further inspection of Fig. 5 reveals that the overall hardness changes are quite 
small, ranging from ~ 0.2 kcal/mol for Si2H6 to ~ 1 kcal/mol for Si2F6. Si2H6, 
Si2Cl6 and Si2Br6 are chemically hardest in the eclipsed conformation, while 
Si2F6 is hardest in the staggered conformation. It is interesting to mention that 
the same trends for the chemical hardness are predicted both by the cation- 
anion energies, at the B3LYP/6-311+H-G** level, and by the HOMO-LUMO 
approximation at MP2/6-311++G** level, as seen in Table 3, where the nu- 
merical values of the hardness in the staggered and eclipsed conformations are 
listed. 



It is noticed that the PMH is verified only for Si2F6, while Si2H6, Si2Cl6 and 
Si2Br6 present hardness profiles obeying the same trend as their energy pro- 
files. According to the PMH, the hardness profile of Si2F6 displays a maximum 
at the stable staggered conformation and a minimum at the unstable eclipsed 
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Molecule 


HOMO-LUMO, B3LYP 


ANION CATION, B3LYP 


HOMO-LUMO, MP2 




0.185 


0.308 


0.52 




-0.995 


-0.856 


-0.90 


SisCle 


0.473 


0.634 




Si2Br6 


0.738 


1.071 


1.39 



Table 3 



Change in chemical hardness Ar] = ije — rjs in kcal/mol. 

conformation. We want to emphasize the complementary behavior of energy 
and hardness: whereas for Si2F6 the almost free internal rotation does not 
allow to distinguish the energetically most favorable a value, the hardness 
profile allows this characterization. In contrast, the hindered rotation in 812X0 
(X=H, CI, Br) yields energetically distinguishable conformations, but they 
cannot be characterized by the hardness profiles. 



6 Reactivity of Silanes 

The reactivity of these systems, induced by the internal rotation, cannot be 
rationalized in terms of the profiles of /i and 77 alone, due to their almost 
constant behavior as a function of a. However, a different perspective of the 
electronic structure and reactivity is provided by the LUMO and HOMO den- 
sities. In Fig. 6 the HOMO of the Si2H6 molecule is shown and we observe very 
similar orbitals to the other 812X6 molecules we have considered. The majority 
of the orbital charge accumulates on the Si-Si bond with some contribution 
on the hydrogens, and with a bond of clear 7r-character. Moreover, there is lit- 
tle difference between the HOMO staggered and eclipsed charge distributions, 
indicating that the effect of the torsional motion on an electrophilic attack is 
negligible. 

The LUMO electronic structure, as illustrated in Fig. 7, is completely dif- 
ferent. For Si2H6, in the lowest energy (staggered) configuration, the charge 
density is delocalized on the sylil groups. Moreover, inspection of Figs. 6 and 
7, shows that the largest overlap between HOMO and LUMO orbitals occurs 
for Si2Fg, which suggests that in the staggered configuration this molecule has 
the strongest hyperconjugative effects (^. On the contrary, in the eclipsed 
conformation the charge is delocalized on the Si-Si bond, with a yr-antibond 
character. This indicates that a nucleophilic attack on Si2H6 may present dif- 
ferent specific mechanisms, as a consequence of the low torsion barrier. In the 
Si2F6 staggered configuration the delocalization process is different: the charge 
is delocalized on the Si-Si bond, with antibond character, and is symmetric 
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Fig. 6. HOMO for Si2H6. Left panel: staggered configuration. Right panel: eclipsed 
configuration 

around the Si-Si bond, but with an asymmetry in the direction of the sylil 
groups. Instead, in the eclipsed conformation, the charge is again delocalized 
on the Si-Si bond with antibond character, but with some preferential charge 
on the sylil side when viewed in a plane with four hydrogen atoms. In addition, 
the delocalization volume is larger in the staggered conformation. 

Finally, we consider the Si2Cl6 molecule (Si2Br6 behaves similarly), which in 
its staggered configuration has the charge localized on the Si atoms, with a 
very clear c-character, but with some asymmetry in the sylil group directions. 
This is similar to the eclipsed configuration, where the charge distributes in 
much the same way, except in that it is completely symmetric around the 
Si-Si bond. The above results suggest that the low, but significant, barriers 
that hinder internal rotation may induce different specific nucleophilic attack 
mechanisms. 



7 Concluding Remcirks 

We have performed DFT calculations on Si2H6, Si2F6, Si2Cl6, and Si2Br6 of the 
evolution of the electronic energies, chemical hardness and chemical potentials 
as a function of torsion angle. For all these molecules at the DFT-B3LYP/6- 
311-1— |-G** level, the staggered conformation is predicted to be the most stable 
one. Moreover, except for Si2F6, it is softer than the eclipsed configuration due 
to a different charge delocahzation at the LUMO orbital. 

Low, but significant, energy barriers hinder internal rotation. For Si2H6 the 
chemical potential and hardness remains quite constant during the torsion 
process, while the other molecules show different degrees of electronic density 
rearrangement as a function of the torsion angle. However, it was not possi- 
ble to characterize precisely the reactivity behavior just on the basis of the 
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Fig. 7. LUMO for the different molecules considered. Left panel: staggered config- 
uration. Right panel: eclipsed configuration. From top to bottom, Si2H6, Si2F6 and 
SiaClg. 

chemical potential and hardness profiles. 

The qualitative analysis of the frontier orbitals shows that for the 812X5 series 
there is little difference between the HOMO staggered and eclipsed charge 
distributions. This indicates that the effect of the torsional motion on an elec- 
trophilic attack is negligible. In contrast, the low but significant barriers that 
hinder internal rotation may induce different nucleophilic attack mechanisms. 
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